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Abstract 

We propose a new mechanism and a related device concept for a robust, magnetic field tunable 
radio-frequency (rf) oscillator using the self oscillation of a magnetic domain wall subject to a uni¬ 
form static magnetic field and a spatially non-uniform vertical dc spin current. The self oscillation 
of the domain wall is created as it translates periodically between two unstable positions, one being 
in the region where both the dc spin current and the magnetic field are present, and the other, 
being where only the magnetic field is present. The vertical dc spin current pushes it away from 
one unstable position while the magnetic field pushes it away from the other. We show that such 
oscillations are stable under noise and can exhibit a quality factor of over 1000. A domain wall 
under dynamic translation, not only being a source for rich physics, is also a promising candidate 
for advancements in nanoelectronics with the actively researched racetrack memory architecture, 
digital and analog switching paradigms as candidate examples. Devising a stable rf oscillator using 
a domain wall is hence another step towards the realization of an all domain wall logic scheme. 
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Self oscillators show a remarkable property of sustaining oscillatory behavior without 
being driven by sources that possess inherent periodicity. In the macroscopic world, a few 
well known examples of self oscillations include the heartbeat, violin string oscillations in 
response to steady bowing^, the Vander Pol oscillator and the infamous collapse of the 
Tacoma Narrows Bridge in 1940^. In the nanoscale too, self oscillations govern the underly¬ 
ing principle of the well studied resonant tunneling diode based oscillator, the spin torque 
oscillator^^® and also, the recently noted phenomenon of nuclear spin induced current os¬ 
cillations in quantum dots®. Our proposal relies on the self oscillations in the translatory 
motion of a magnetic domain wall. 

The interest in dynamics of magnetic domain walls has been active for decades^^^^, re¬ 
cently intensihed by the discovery of current-driven domain wall motion^^^^® and its related 
applications in nanoelectronics®’^®^^®. Field driven oscillations have also been observed and 
studied^^’^® for a long time. However, these oscillations are accompanied by a drift which 
makes them unusable as a device. On a different note, vertical injection of uniform spin 
current is proposed as a means for high domain wall velocities^®^^^. Here, we propose stable 
oscillations caused by a constant magnetic held whose drift is canceled by a vertically in¬ 
jected non-uniform spin current as depicted in the schematic in Fig. 1(a), thereby resulting 
in a stable periodic motion. 

Normally, domain wall motion under vertical spin currents is caused by field like torque 
which is typically smaller than the Slonczewski like torque^'^^'^^ usually responsible for the 
switching of the free layer of magnetic tunnel junctions (MTJ)^®. However, the presence of 
a magnetic held allows an efficient transfer of the Slonczewski like torque without any need 
of a held like term, as shown later. We note that current-driven domain wall oscillations 
with a drift have been observed^^ while stable oscillations on a pinned domain wall were 
proposed^®. Our proposal dihers from the latter in the control of oscillations via an external 
held instead of a pinning potential, both of which under rigid domain wall approximation 
work to rotate the domain wall. 

We hrst provide a theoretical analysis using a rigid domain wall model to hnd an ap¬ 
proximate waveform for the oscillations. We show that the frequency is twice the resonant 
frequency of a magnet in a magnetic held; while the amplitude is approximately a linear 
function of the ratio of the hard axis anisotropy and the magnetic held. The oscillatory 
part of the waveform is independent of the input spin current to a very good degree of accu- 






FIG. 1. Working principle and device design, (a) Schematic depicting the principal idea of the 
proposal which comprises a domain wall magnet with the incident non-uniform spin current and 
a uniform magnetic field. The domain wall is shown at an arbitrary instantaneous position. The 
non-uniform spin current is incident on a region spanning one-half of the magnet, (b) A 3D 
schematic of the proposed device design. A strip of Tantalum localized in the desired region of the 
domain wall magnet, connected to a current source generates the desired spin current. The applied 
magnetic field along the —axis is not shown here for clarity, (c) A circuit diagram representing 
the device schematic in (b). ri is the resistance of the part of Tantalum strip between the constant 
current source, ‘T, and MTJ, and similarly r 2 - (d) Zoomed-in motion of middle spin marked in (a), 
depicting the rotation it undergoes which gets converted into motion via its hard axis anisotropy. 
The orange colored region is the portion where spin current is non-zero. 
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racy; and hence can be of great technological advantage for accurate oscillatory waveforms. 
We then also numerically simulate the micromagnetic Landau-Lifshitz-Gilbert equation (see 
supplementary material) as a verihcation of the oscillations and hnd the frequency to be in 
good agreement with the analytic result. Lastly, we analyze the effect of thermal noise on 
the oscillator and hnd that it is robust at room temperature by numerically calculating its 
Q-factor. 

The working principle of the proposed oscillator is depicted in Fig. 1(a). The vertically 
incident spin current is spatially conhned to the left half of the domain wall, while the uniform 
magnetic held along the — 2 : direction exists throughout the region. Both the regions, 2 : ^ 0 
and 2 : 0 by themselves are unstable for the domain wall (unless the spin current is too low) 

and hence it is restricted to be in some region around z = A possible device realization 
of the proposal, in which the domain wall self-oscillations may be ehectively translated into 
an alternating current oscillations is depicted in Fig. 1(b). A dc-current source is connected 
to the bottom strip of Tantalum, which is used to inject a spin current via the giant spin 
hall ehect (SHE)^®^^^. The strip of Tantalum conveniently acts as a source of spin current 
that injects a constant spin current in a locally conhned region marked in the schematic in 
Fig. 1(a). An MTJ like structure is then used to sense the position of the domain wall via a 
measurement of the change in resistance^'^. A small current is applied and the corresponding 
voltage across the MTJ is measured. An equivalent circuit for the entire set up is shown in 
Fig. 1(c) depicting the measurement in a more explicit way. We assume that the current 
used for the measurement is small enough such that it doesn’t have any additional ehect on 
the domain wall dynamics. 

The analysis to follow will be based on Fig. 1(a) which captures the essence of our 
proposal. The region of non-zero spin current “pushes” the domain wall towards the region 
of zero spin current via spin transfer torque (STT). However, the domain wall cannot keep 
moving away from the spin current as the magnetic held will push it back via the dissipative 
Gilbert term®. As shown in Fig. 1(d), when the domain wall enters the spin current region, 
it is rehected back with a diherent azimuthal angle accumulated because of the magnetic 
held. The hard axis anisotropy then keeps the domain wall moving until the magnetic held 
rotates the domain wall again to cause reverse motion. Hence, the magnetic held causes a 
perpetual rotation, while the hard axis anisotropy converts the rotation into a translation 
of the domain wall. The spin current then acts as the energy input which negates the 
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dissipative effect present in a purely field driven motion, hence stopping the drift observed 
in the latter case. For a typical case of a low dissipation constant, the domain wall will need 
only a small amount of push from the spin current and hence the oscillations will be almost 
independent of it. The average location of the domain wall though will be dependent on the 
spin current density. Our analysis is based on the magnetization dynamics of the domain 
wall described by the Landau-Lifshitz-Gilbert equation augmented by the Slonczewski spin 
torque term^^ given by 

— = - 7 m X (iJ„. + H.,) + X (^— j - X (J [z) X m), (1) 

where m{z,t) is the magnetization unit vector; 7 (> 0) is the gyromagnetic ratio; Mg is 

the saturation magnetization of the magnet; /tq is the permeability of free space; d is the 

thickness of the sample (Fig. 1(b)); a is the Gilbert dissipation constant; J^{z) is the vertical 

spin current density loss which is assumed to be dependent only on ; 2 ; iText is the externally 
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applied held; iTeff = —~ H±my + is the exchange energy constant; H± 

is the hard axis anisotropy; is the easy axis anisotropy. The above equation can be derived 
from the Lagrangian along with the generalized forces (see supplementary information) given 
by the following expressions. 
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where w is the width of the magnet (see Fig. 1(b)). We consider the rigid domain wall 


ansatz^^, (l){z,t) = and 9{z, t) = 2tan ^ exp 
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We consider the case of a spin 


domain wall, A(^/>), is given by /zj , rj ■ 2 

fioMg [H\\ + H± sm^ Ip) 

current polarized only along the G’ axis with its expression being J® = —Jp6{—z)z where 6 
is the Heaviside function. Additionally we have a uniform magnetic held along —direction, 
Hext = —Hz. Finally, we get the equation of motion as. 
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where the dimensionless time is dehned as r = 2'yHt and we dehne constants G± = —— and 

2H 

flj^ 

Gi =- I . The frequency can directly be deduced from eqs, (4) and (5), as follows. 

2ejj,QMsHd 

Let there exist an oscillatory solution of Z and A with a common dimensionless time (r) 
period, say T. Note that under such conditions, all the terms in Eq (4) are periodic with 
period T except possibly the term containing sin cos V’- For this term to be periodic with 
T, tfj has to change by an integer multiple of tt after the time period. Considering Eq (5), 
it can be seen that ip has an oscillatory term along with a drift of “speed” —1/2. Using 
the known value of “speed” of ip, we can conclude that any oscillations have to exist with a 
dimensionless time period of 2 n 7 r where n is any natural number. We restrict our attention 
to T = 27r which is what we observe in simulation. Coming back to real time, t, we conclude 
that rigid domain wall approximation restricts the angular frequency of oscillation to be 
2'yH (see Eq (7)). Note that with this time period, we indeed have an oscillation in A with 
the same dimensionless time period. 

If Gi = 0, the equation for ip can be solved exactly® and Z can be integrated using 
Eq (5). With Gi 7 ^ 0, there does not seem to be an analytic solution. However, using the 
intuition that spin current is a small perturbation which mainly acts to negate the effect of 
dissipation, we can use the solution for held driven motion® to approximate the oscillations 


as, 


Z(r) = Zc(Gi,Gx) + ^ 

a 


_ T\ IT — TT 

tan“^ ( q;G_l + tan - ) + 27rCeil —-— 
2 / \ 27r 


( 6 ) 


where Ceil is the ceiling function; Zq is a constant dependent on spin current and helds 
while the oscillatory part is dependent only on G_l; (A) is the average value of A which 
can be approximated as \{ip = 7 r/ 4 ). We can now use Eq ( 6 ) to analyze the waveform 
explicitly to calculate various useful observables. First, we use it to hnd the amplitude of 

We also note that for q;Gj_ 1, this expression reduces 


• 1 . . 2(A) aG 

oscillation as -tan 


a 2 

approximately to Eq ( 8 ). The fact that hard axis anisotropy is responsible for converting 
rotation into motion is backed up by the expression of the amplitude of the oscillation being 
approximately proportional to the ratio of the hard axis anisotropy and the magnetic held 
(see eq 8 ). To calculate the threshold of spin current, we take the time average of Eq (4). 

/ Q'i ' 

Then using the waveform (see Eq ( 6 )) in Eq (5), we get the inequality, Gj > a f 1 H—^ 
This is in accordance with the intuition that the spin current mainly acts to negate the 
ehect of dissipation present in the held driven motion of a domain wall. The results after 
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FIG. 2. Oscillatory waveform and region of occurence (a) Comparison of the oscillatory part of 
the numerical and the analytic solution of domain wall position vs time, under an external applied 
field of lOfcA/m. The spin current density for the numerical solution was taken to be 0.96GA/m^. 
(b) The values of spin current and magnetic field depicting the region in which the oscillations 
happen. 

dismantling the notation are summarised in Eq (7), Eq (8) and Eq (9). 
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where \eq is the equilibrium value of A and the factor in the square root (Eq (8)) arises 
due to the variation in the width of the domain wall under oscillatory motion. It will be 
absent if the variation in width is neglected or in other words, <C H\\. 

We now demonstrate the simulated results of the domain wall motion using the rigid wall 
approximation discussed above. The waveform derived in Eq. (6) is not an exact solution of 
the equations for the rigid domain wall, but matches fairly well with the numerics as shown 
in Fig. 2(a). The regime of operation of the device as shown in Fig. 2(b) demonstrates that 
we need a minimum magnitude of the spin current to compete against the magnetic held 
and hence result in the oscillations. This can be understood by analyzing the motion of the 
domain wall when it starts from deep inside either region, i.e., the region deep inside the 
region of zero or non-zero spin current. 

In Fig. 3(a), we demonstrate the simulated motion of a rigid domain wall starting from a 
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FIG. 3. Current push and Field push, (a) Simulation results for domain wall position, Z when 
started from a point inside the region of non-zero spin current, (b) Similar plot as (a) except for 
the initial position of domain wall being inside the region of zero spin current. 


point < 0. As shown in the hgure, the domain wall will be pushed away until the “force” 
of the spin current is small enough to be compensated by the drift caused by the magnetic 
held. An opposite scenario is shown in Fig. 3(b), where the domain wall starting deep inside 
the region of zero spin current {z > 0) will have a held driven drift until it encounters the 
region of non-zero spin current. In both the scenarios, the spin current magnitude should be 
large enough to push back the domain wall or the latter will continue to move indehnitely 
against the held. 

To verify the results, we have also performed micromagnetic simulations, details of which 
are included in the supplementary material. We consider a 3nm thick magnetic him with a 
cross-section of 800 x lOOnm. We assume the magnet parameters. Mg = 8 x 10^ A/m and 
Aex = ISpJ/m; where Mg is the saturation magnetization and A^x is the exchange energy. 
We assume a crystalline anisotropy, with its contribution to energy density as —Karriy where 
Ka = 0.35MJ/m^, in the direction of its thickness which works toward reducing the hard axis 
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FIG. 4. Noise Analysis. Power Spectral Density of the domain wall position for three values of 
the applied field and with the spin current density hxed at 0.96GA/m^. It shows peaks at the 
expected values of the oscillation frequency with small width due to noise. 


anisotropy caused by the dipolar interaction. We apply a magnetic field of 8.75 kA/m and 
a spin current density of 20 GA/m^. We have simulated it for 40ns to hnd that oscillations 
occur at a frequency close to 0.56GHz. This value is close to the one derived using the 
theoretical analysis as written in Eq (7). 

Noise Analysis: Finally, we verify the stability of the oscillations by adding noise to the 
rigid domain wall equations and numerically calculating its quality factor. For simplicity, we 
assume H± which corresponds to a domain wall with constant width. We introduce 

two uncorrelated white noise sources, Nz and N^, both of which satisfy {N(t)) = 0 and 
{N{t)N(t')) = 2akBT6{t — t'). The noise can then be added as^^. 


IdZ 
A dr 
d^jj 
dr 


^ , , di/j Gi 1 

Gxsm'ipcos'ip + a— + —-+ 

dr 2 1 + exp (^) 

1 IdZ W 

— — (y - 1 - ^ 

2 A dn Gn 



( 10 ) 

( 11 ) 


where Gn = A’ynQMgH'^wdX has been chosen such that Fokker-Planck equation correspond¬ 
ing to Langevin equations Eq (10) and Eq (11) admits the Boltzmann distribution in steady 
state. We simulated the above equations at a room temperature of 30077 for 40/rs for various 
values of spin current and magnetic held. The power spectral density (PSD) of Z for three 
values of the applied magnetic held is plotted in Fig. 4. From the simulated spectrum we 
hnd that the quality factor of the oscillator is ~ 550, ~ 1100 and ~ 1400 respectively for 
the applied helds of 8kA/m, lOkA/m and 12kA/m respectively. 

In conclusion, we have proposed a new set up for an oscillator based on the self oscillations 
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of a magnetic domain wall. We found that under rigid domain wall approximation, the 
oscillatory part of waveform is almost independent of input spin current and the frequency 
of oscillations is solely governed by the external magnetic held. We also demonstrated a 
high quality factor giving evidence for the stability of the oscillations. We envision that 
the simple set up proposed, namely a domain wall subject to a non-uniform vertical spin 
current will also open up many possibilities for simultaneous write and read out along with 
the possibility of an all domain wall logic scheme. 
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Abstract 

Here, we provide details about the rigid wall approximation along with the details of the micro- 
magnetic simulations performed in the paper. 
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DERIVATION OF LAGRANGIAN 


The starting point of magnetization dynamics is the Landau-Lifshitz-Gilbert (LLG) equa¬ 
tion augmented with Slonczweski spin torque^ term 
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dt 


= — 7 m X ifeff + Oim X 


dm 

dt 


'jh 


efioMsd 


-m X (J^(z) X m) 


( 1 ) 


where J^{z) is the incident spin current density which is assumed to be dependent only on 
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from the Lagrangian, eq (2), along with the generalized forces derived from the expression 
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where fioMgE = E. In this section, we derive this equivalence. We take {mz{z),(l){z)} as 
coordinates. Hence, we can write. 


5m = sin 0(50$ -|- 59& 


(4) 


where cos 9 = niz, $ = ( — sin 0, cos 0, 0) and © = (cos 9 cos 0, cos 9 sin 0, — sin 9). This gives 
m.5m = sin^ 00(50 -|- 959 which translates to, 

m.5m = sin^ 9^5(1) + (^ 5 ^ 

1 — mj. 

One can notice that {m, ©,$} forms an orthonormal right handed system. They are 
just the unit vectors in polar coordinates with m acting like the radius vector. Using 
Euler Lagrange equations by varying the action with respect to the coordinates, we get two 
equations of motion. 
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The above equations are equivalent to the Landau-Lifshitz Gilbert (LLG) equation. To 
show this, we explicitly reduce down the terms in eq (6) to relevant terms in eq (1). 
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From the above set of equalities, we can readily write eq (6) as, 

^ f dm\ 'yh . , -rsf \ 

— = -Um X + a[rn y — 'j- (m x (J (x) x m)). (15) 

which is just the z-component of eq (1). Hence, one of the equations we get from Lagrangian 
is compatible with LLG. For the next equation, eq (7), we proceed from eq (1) and show 
that we can derive eq (7), 
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Hence, we show that both eq (6) and eq (7) are satished if the original LLG equation 
eq (1) is. As the number of independent equations in LLG is 2, we can say that the former 
two equations are equivalent to LLG. 
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DERIVATION OF THE RIGID DOMAIN WALL APPROXIMATION 


We concluded the last section by deriving a Lagrangian for the magnetization dynamics. 
In this section, we use rigid domain wall approximation to derive the equation of motion for 
the collective coordinates. This amounts to putting the following ansatz in the Lagrangian, 


0(x,t) = V’(t) 

(17) 

0{z,t) = 2 tan ^ (e j 

(18) 


After we put this ansatz in the Lagrangian eq (2), we can carry out the integration with 
respect to the 2 ; co-ordinate. The energy functional (divided by fioMg) is, 

E = {d,mf + - ^ml - m.ifext (19) 

Then, we will be left with a Lagrangian with Z, and A as coordinates. It is this hnal 
Lagrangian that will give us the equation of motion via the Euler-Lagrange equations. The 
generalized force due to STT, however, will have to be changed as well. By the prescription 
of work-energy, we have 

J dz {r X m) .5m = + Fz5Z + Fx5\ (20) 

where are functions of {Z,'^,A} which we are going to hnd. For convenience 

in notation, we resolve as. 


J®(x) = Ji{z)m{z) + J 2 {z)m'{z) + Js{z)m{z) x m'{z) (21) 


where m'{z) denotes the derivative of m{z) w.r.t x. The above expansion is valid only if m' 
is non-zero everywhere. This is true in the particular ansatz we are considering. We write 
down the relevant expressions now. 


m'. (m X 5m) 

m' .5 m 


- sin^ e 


5il) 

T 


- sin^ 0 


(5Z 5\{z-Z)\ 

f V ^ v ) 


( 22 ) 

(23) 


Now, using the above equations, the term under consideration is (refer (21) for dehnition 
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of Jj), 


. {m X 5m) = J 2 m'.{m x 5m) + J^m'.5m 
= —J29'sm.95'il) + J'i9'59 
= -sin\9{z)) 


J2{z) 


5i){t) 


m 


dzJ\ (m X 5m) = F^5ij + Fz5Z + Fx5X 

f dz 4 

— — 

Fz=- 
Fx = - 


[exp (^) + exp (- 


dz 4 


[exp + exp (- 


dz 4(x — Z) 



[exp (^) + exp (-^) 



(24) 

5Z 5X{z-Z)V 
A 2 + A3 )_ 

(25) 

(26) 


(27) 

2-Hz) 

(28) 


(29) 


(30) 


So far, it was generally applicable for any distribution of spin currents. For the case under 
consideration, we focus on a specific spin current distribution such that is a constant, say 
J|, for X < 0 and zero for x > 0. and J'^ are zero everywhere. Then, we get Fz = Fx = 0 
with given by. 


1 + exp 


(31) 


Now, we calculate the integrals with respect to the x co-ordinate for the terms in the 
main Lagrangian. We list down the equality among integrals. 


dz- 


-m. 


7 


-2Z . 

-^ 

7 


A 9 A 

J ^ex / r., \2 

dz — —{^x'>n) = ——— 

fioMg fioMgX 

j dz^^{m.y)‘^ = H±Xsm‘^^jJ 

[ dz - -{m.z)'^ = Htdl — X) 


'a . ^ 2aX (Z5Z 
dz—m.5m = 


7 


7 \ A2 


+ i>5ip + 


\ 


crTT^ X5X 
67 A 


dz (—m.Hp^t) = —2Hp^tZ 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 
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The final Lagrangian along with work is given by, 
^-‘IZ ■ 2 A. 


L = noMsiud 


5W = —fioMswd 


7 




HoMgX 


i7_i_A sin^ -i/j + (/ — A) + ttA {H^ cos il) + Hy sin ip) + 2HzZ 


2a\ Z.5Z 


7 


A 2 


+ ip.6ip + 


ttTT^ \.5\ 
67 A 


+ Fz6Z + F^6iP + FxSX 


(38) 

(39) 


The Euler-Lagrange equations then lead to, 

; TT 

-Ip + qiiext = ^ 


-— 'yH± sin Ip cos Ip = aip 

A 




enoMgd 1 + exp (^) 


2 A 


— H± sin^ Ip — H\t = 


aTT^ A 


(40) 

(41) 

(42) 


Note that we have applied a field in —direction in the paper and hence in the notation of 
the paper, ifext = —H- Similarly, the sign of Jp here is opposite of that in the paper. 

So far, we have treated all coordinates {Z, ip, X} on equal footing. However, A is different 
from the other two in the sense that it changes the shape of Domain Wall instead of just 
coherent translation (Z) and rotation [ip). Hence, under the validity of the rigid domain 
wall, the variation in A should be small. This can be verified through numerical simulations 
as well. Using this approximation as a guide, we neglect the time derivative term (which 
is also proportional to a, thereby helping the approximation) in eq (42). Hence, we arrive 
at the equation of motion in {Z,ip} with A as a function of ip instead of an independent 

variable, given by X{ip) = 


2 A 


^oMs{Fl\\ + TT^sim V’) 


NOISE ANALYSIS 

At T = 0, rigid domain wall approximation correspond to eqs (40) and (41) with qualifi¬ 
cations about A from the subsequent paragraph. At non-zero temperature, we have to add 
terms pertaining to noise in the aformentioned equations. In this section, we derive these 
terms and discuss the simulations of the stochastic differential equations. For analytical 
simplicity, we only consider the case when H\\, such that A can be taken to be a 

constant. Also, we consider the case of equilibrium without an external magnetic field or 
spin current and then assume that the noise sources remain the same in non-equilibrium. 
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Under such circumstances, we can rewrite the hnal set of equations more generally by the 
following dehnitions, 


P = 


X = 


12iiQMsWd 

7 

! 2ij,oMswd 

7 




S = noMsH±wd\ sin^ V’ 


(43) 

(44) 

(45) 


Note that with A as a constant, the only term remaining in energy is the hard axis anisotropy. 
After some calculations, we can reduce eqs (40) and (41) in terms of x and p as. 


X = 


P = 


1 + a2 
1 


dS ,dS\ 

o n ) 9x^x 

op OX J 


_ (d^ ad^\ 

1 + \dx^ \ dp) ^ ^ 


(46) 

(47) 


where we have added two uncorrelated white source, and Ap; and Pp are assumed 

2q(/c 99 

to be independent of Z and ip. We assume {N(t)N(t')) = - -6(t — t') for both and 

1 + 

Np. These assumptions are justihed by deriving the expressions of and gp from the 
Fokker-Planck equations corresponding to the above Langevin equations. They are given 


— = toWpl - — (pWpl + — rp(2)pi p ^ rp(2)pi 
dt dx^ ^ J dp^ P J 9x2 L - Qp2V P J 

where the drift and diffusion coefficients are given by, 

£)(i) = —1— - aX—'] D<'> =-i— (— + 

* l + aH9p OxJ <• l + a'‘\dx Xdp) 

Df! = akaTgl = aksTc/^ 


(48) 


(49) 

(50) 


After a series of elementary manipulations, we can get the expressions for g^ = \ and 
gp = 1/A by demanding that P = C eyip{—£/ksT) satishes the Fokker-Planck equations 
in equilibrium, where C is some constant. Now we revert back to the original coordinates, 
{Z,We dehne two noise terms, Nz = N^. — aNp and W/; = oiNg. -f- Np to simplify the 
hnal expression as. 


A 

ip 


■yH± sin ip cos ip + aip -\- 


jhJp 


TT 

iHext I" 


7 


2p,QMsWd\ 


epoMgd 1 -|- exp(2Z/A) 


+ 


7 


2fj,oMswd\ 


-N. 


(51) 

(52) 


7 



where we have added the external magnetic held and spin current back. Note that both Nz 
and are uncorrelated white noise sources with {N(t)N(t')) = 2akBT6{t — t'). 

We then numerically simulated eqs (51) and (52) using Stratonovich calculus and second 
order Heunn scheme. We wrote the code in MATLAB by discretising time in units of 
Ips and assuming the parameters as Mg = 800A;A/m, Aex = ISpJ/m, = 32.9kA/m, 
H± = 16.QkA/m, 7 = 2.21 x 10^m/(A — s) and a = 0.007. To verify the code under 
equilibrium, we added a pinning potential of the form OAkpi^Z"^, with kpi^ = 15.12/iJ/m^, 
and assumed a low temperature of 5077 such that sin^'^ ^ Under such circumstances, 
the energy (refer eq (45)) becomes analogous to a particle in a harmonic potential and we 
verihed that the equipartition theorem holds by hnding the average value of and 
All averages as well as any other statistical values were calculated after removing the initial 
motion for 40ns to allow the system to achieve time invariance. 

We then applied the external inputs and assumed a room temperature of 30077. We 
simulated the motion for 40/is and found the spectral function (again after removing the 
initial 40ns of motion). The simulation was run for various values of 77ext and and we 
extracted the Q-factor in all the cases, which turned out to be between 500 and 1500. The 
simulation result for three cases is shown in the paper (see Fig 4 of the paper). 


MICROMAGNETIC SIMULATIONS 

The theoretical analysis in the main paper was carried out using rigid domain wall ap¬ 
proximation imposed on a 1 dimensional LLG with effective parameters for anisotropy. For 
a more quantitative treatment, we must take into account the dipolar interaction precisely 
instead of masking it with a local effective held as in eq (19). This, in general, might be a 
difficult problem to solve analytically. Hence, we simulate a magnetic domain wall under 
the conditions given in the paper to verify if the oscillations are happening and compare the 
motion with the one derived from theoretical analysis. 

We consider a magnetic thin him of thickness 3nm (along y-axis), width lOOnm (along 
x-axis) and length 800nm (along x-axis). Consider the LLG equation given in eq (1) with 
the ehective magnetic held given by a sum of exchange held, crystalline anisotropy, held due 



to other dipoles and an external magnetic field. 


2Aex rj2 , 

——-— 


fj,oMs 


fj,oM, 


■rriyy + i?ext + i?dip 


where we assume crystalline anisotropy, Ka, to reduce the hard axis anisotropy due to 
demagnetisation field (this is required for rigid domain wall approximation to work); i?dip 
is the effective field due to the dipolar interaction between magnetizations. i?dip has been 
calculated as discussed by Miltat et al.^ 

We use finite difference method and fourth order Runge Kutta with a fixed time step of 
0.2ps to simulate a tail-to-tail domain wall. As the thickness of the film is small compared to 
typical exchange length, we perform a 2-dimensional simulation with discretization in blocks 
of 4nm X 4nm. The values of parameters are taken to be Agx = Mg = 800kA/m, 

Ka = 0.35MJ/m^ and a = 0.007. To achieve the domain wall configuration, we simulate 
an infinite strip by adding additional magnetic charges at the two edges along the length^. 
These magnetic charges are due to the remaning left and right parts of the infinite strip, 
which are assumed to have a constant magnetization pointing along ±x-axis. Under such 
conditions, we let the system equilibriate for a sufficient amount of time. After equilibrium, 
we observe that my is zero in the bulk and develop small non-zero values at the edges 
along the width. We average {mx,my,mz} along the strip width (x-axis) to get the plots 
given in Fig 1. We have also shown the theoretical waveform derived from rigid domain 
wall approximation, mz{z) = tanh ^ with the domain wall position, Z and width, 

A chosen to fit the micromagnetic result. The system equilibriates ai Z = 64.2nm and 
A = 38.9nm. We note that the initial location of the domain wall, Z, is dependent on the 
initial conditions of the simulations while the width. A, is independent of it. 

We then apply an external magnetic field of 8.75kA/m and vary the vertical spin current 
density to see the motion under these stimuli. The spin current applied is restricted to 
one half of the magnet as discussed in the paper. Fig 2 shows the motion under various 
values of spin current density, all in GA/m^. The domain wall position has been defined 
as {—1/2) * {mz)x where the average is taken along x-direction. Note that we could as well 
define Z to be the point where {mz)x is zero. However, this distinction doesn’t quantitatively 
affect the dynamics. At J = 0, we get back the usual oscillatory and drifting motion of the 
domain wall position. For non-zero small J, we still have a drifting domain wall but with a 
smaller velocity. For higher J, we get pure oscillations with the average domain wall position 
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FIG. 1. The equilibrium magnetization along the length of the magnet. Here {mx, my, m^} denotes 
the average of the quantities taken along the width, x-axis. The dots (‘RDW’) refers to the fitted 
waveform of in rigid domain wall approximation. 


100 



FIG. 2. The motion of the domain wall under various values of spin current density with a magnetic 
field of 8.75kA/m. The spin current is restricted to be in one half of the magnet as discussed in 
the paper. 

dependent on the spin cnrrent as discnssed in the paper. 

Now, we focns on the case of J = 20GA/m‘^ where oscillations are stable. The oscilla¬ 
tion freqnency is 0.56GHz which agrees with the one derived from the rigid domain wall 
approximation in the paper, 0.61GHz. As discussed in the paper, we have both translatory 
oscillations as governed by Z and rotation of the wall as governed by V’- To show the rotation 
we plot the time evolution of {my)x and {mx)x at the instantaneous domain wall position in 
Fig 3. We can see that they are oscillating with a phase difference of close to 7r/2. In rigid 
domain wall approximation, the oscillations of 'll) occur with half the frequency of Z. This 
can be explicitly noted in micromagnetic simulations as well by plotting {my)x at z = Z{t) 
as a function of time and comparing it to Z{t) itself as shown in Fig 4. From eq (41), we 
can see if the domain wall is sufficiently outside the region of non-zero spin current, Z ^ 
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FIG. 3. Depicting rotation in domain wall. Here niy means the average value of my at the 
instantaneous position of the domain wall. The averaging is done across the width, is computed 
similarly. 



FIG. 4. Illustrating the phase relation between the domain wall position and its rotation. Here 
again, my denotes the average value across the width at the instantaneous position of the domain 
wall. 

the velocity of the domain wall is zero at V’ = 0, ±7r/2,7r. Even thongh this condition is 
never achieved in the shown oscillations, it can be approximately noticed from the plot at 
the maxima of the domain wall position (where Z is still less than A). 

Fig 5 is a plot of {mz)x at the fonr points marked as {A, B, C, D} in Fig 4. These 
waveforms can be htted with rigid domain wall approximation with the domain wall position 
and width as htting parameters. Using this, we can verify that at least for the averaged 
value of TTiz, the rigid domain wall approximation is valid throughout the oscillations. This 
also illustrates the oscillations of the width of the domain wall. 

To compare the micromagnetic simulations with theoretical analysis, we need to dehne 
H± and through micromagnetics. The dehnition is given in terms of mean held of a 
mono domain magnet. We extract as the mean held when all the magnetic moments are 
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FIG. 5. The plot of ruz averaged across the width for four different instances in the oscillations. 
This depicts the variation in width during different phases of the oscillations. 

pointing along x axis. Hy and Hz are defined analogously. We then define H\\ to be — Hz 
and if_L as Hy — H^- We note however that dehning a particular H±^ and H\\ can only give 
an estimate of the observables and cannot be used for quantitative correctness^. We get 
the values as H\\ k, 32.9kA/m and H± ^ 16.6kA/m. Using these as effective parameters 
in the LLG equation, we calculate the amplitude of oscillations through the rigid domain 
wall approximation. We hnd that the amplitude is of the same order as the micromagnetic 
prediction. The amplitude in micromagnetic simulations is approximately twice of the one 
derived from theoretical analysis in the paper. The factor of 2 could be because of a specihc 
dehnition of H_i and H\\ and hence cannot be predicted easily. This is an evidence of the 
validity of the theoretical analysis carried out in the main paper. 
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